High Energy Neutrinos from Gamma-Ray Bursts with Precursor Supernovae 
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The high energy neutrino signature from proton-proton and photo-meson interactions in a super- 
nova remnant shell ejected prior to a gamma-ray burst provides a test for the precursor supernova, 
or supranova, model of gamma-ray bursts. Protons in the supernova remnant shell, and photons 
entrapped from a supernova explosion or a pulsar wind from a fast-rotating neutron star remnant 
provide ample targets for protons escaping the internal shocks of the gamma-ray burst to inter- 
act and produce high energy neutrinos. We calculate the expected neutrino fluxes, which can be 
detected by current and future experiments. 

PACS numbers: 96.40.Tv,98.70.Rz,98.70.Sa 



Gamma-ray bursts (GRBs) are thought to be possi- 
ble sources of high energy neutrinos. In the currently 
favored models the 7-ray emission is attributed to radi- 
ation from shock-accelerated electrons in the relativistic 
fireball outflow or jet from a cataclysmic stellar event. 
The latter may be connected to compact stellar merg- 
ers, or the core collapse of a massive stellar progenitor 
(collapsar), which could also involve a core-collapse su- 
pernova (SN) 0]. Together with electrons, the shocks are 
expected to accelerate protons as well, and high energy 
neutrinos are thought to be produced dominantly by pry 
interactions of the protons with synchrotron or inverse 
Compton scattered photons 0, • Recent reports of de- 
tection of X-ray lines from several GRB afterglows have 
been interpreted Q , although not unambiguously, as pro- 
viding support for a version of the collapsar model in 
which a SN explosion occurs weeks before the GRB (the 
"supranova" model, [f|). In the supranova scenario the 
supernova remnant (SNR) shell provides nucleon targets 
for pp interactions with protons accelerated in the MHD 
wind of a pre-GRB pulsar 0, leading to a 10 TeV neu- 
trino precursor to the GRB (other nucleon targets from 
a stellar companion disruption leading to ir° decay GeV 
7-rays were discussed by 0). The SNR also provides 
additional target photons for pj interactions [6j with in- 
ternal shock-accelerated protons, resulting in ~ 10 16 eV 
neutrinos. In this Letter we investigate some important 
unexplored aspects of the GRB-SNR interaction, namely 
pp and pry interactions involving GRB protons in the 
shocked shell, which have significant consequences in the 
assessment of the neutrino signatures from these objects. 

SNR target nucleons and photons. — The typical mass 



ejected in a SN is M S] 



r M m = 2 x 10 33 r 



g. For 



a nominal sub-relativistic shell speed v = 10 cm/s, the 
typical distance reached is R snl ~ 10 14 fgi ( j cm in td days. 
We assume that the gas is roughly isotropically ejected in 
a shell of width 6 = AR/R = 0.1<S_i and average column 
density £ « 1.6 x 10 4 m s 
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tor stellar envelope of dimension i?„ = 10 12 i?i2 cm 
heats it to a temperature T Q ~ 1 E\J^ 8l\ keV. 
The mean photon energy in the SNR shell is e 7iSn ~ 
2QEl[ i 5 1 _!^R\I^VQ 1 t~ 1 eV. The photon column density 
in the SNR is £. 



7,sn 
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cm . The P7 optical depth and threshold proton en- 
ergy at the A-resonance are 
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10 r E^ /4 SZl /3 Ri 2 1/4 V9t d GeV. (1) 



The collapse of the Fe core can lead initially to a 
fast-rotating pulsar of rotational energy E rot ~ 10 53 



and spin-down time ~ 3 x 10 s 



erg 

|, emitting an 
MHD wind of luminosity L m ~ 10 46 L m 46 erg/s. Be- 
sides ejecting the SNR envelope, the SN explosion may 
leave behind parts of the He core, which take longer in 
falling back to make a black hole leading to the GRB 
(months in the supranova scenario). If the fall-back is 
anisotropic and a channel forms in the He core (e.g. 
along the rotation axis), the relativistic MHD wind may 
flow out and impact the SNR shell further out. The 
MHD wind is highly relativistic, and is decelerated to 
a subrelativistic velocity in a reverse shock, driving a 
forward shock into the shell The shock velocity v s 
in the observer frame satisfies, from pressure balance in 
the shock frame, v s - v ~ [L ra (lQ- 1 R/c)/M sm ] 1 / 2 ~ 
lO-^tVMsnr] 1 / 2 - 10~ 2 c for typical E{t) = tL m 
and M snr . Thus the shock will propagate through a 
significant fraction of the shell width, and we approxi- 
mate the situation by taking the shock to cross the en- 
tire shell, producing photons distributed throughout the 
shell. Given the high Thomson optical depth in the SNR 
shell, these photons thermalize to a maximum energy 
£7,111 ~ 0.1 Lm4 6 6_iVg 3 ^ 4 t d keV, assuming that all 
the MHD wind energy goes into photons. The pj optical 
depth and threshold proton energy at A production, are 



The deposition of 10 E51 ergs in a SN progeni- 
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E p .t Km ~ 2 x 10 6 L^<C 1 1/ V 3/ V /a GeV. (2) 

Additional target photons may arise from parts of the 
disrupted He core, no longer in hydrostatic equilibrium 
and moving outwards inside the shell. However, even 
if they have super-Eddington luminosities L^Hc ~ 10 44 
erg/s as inferred in some SN, these photons have an op- 
tical depth t p1 ~ bL^Vg 2 ^ 1 for protons with E p > 
3 x 10 8 td GeV, negligible compared to Eqs. and (0). 

The non-relativistic MHD forward shock moving into 
the SNR is likely to be collisionlcss (as also generally 
assumed for non-relativistic GRB internal shocks [lflj). 
The transition is dominated by the process giving the 
narrower shock width, i.e. the shortest timescales for 
changing particle momenta. A collisionlcss shock thick- 
ness is ~ c/(/3u) p ), where f3 is shock velocity, u> p is plasma 
ion frequency. For our parameters ~ 1 M Q shell at 
- 10 14 cm radius, u p = 3 x 10 10 /sec and (3 = 0.01, 
this shock width is 100 cm. The thickness of a radiation 
shock is 1/(3, where I is the photon mean free path, giving 
a thickness ~ 10 11 cm. The thickness for pp collisions is 
even larger due to a smaller cross section than Thomson. 
This strongly suggests that the shock is collisionless. 

The average magnetic field inside the SNR shell is 

B ~ 10 4 d 7 VV r 2 ^ V 7 V /2 G, (3) 

where £g is the equipartition value in respect to the 
proton kinetic energy. Equating the acceleration time 
t a ~ 2irAmc^/{f3 2 eB) [where = (v s - v)/c ~ 10~ 2 /?_ 2 
is the relative shock speed in the shell, and A ~ lCL4i] 
to the synchrotron loss time t sy — 3m 3 c 5 /(2e 4 B 2 j), 
the maximum accelerated particle Lorentz factor is 
7mx = 2 x 10 10 A" 1 ' 2 '(m/m p )B-^ 2 . For the field of 
Eq. ©, the maximum proton Lorentz factor is j Pt mx ~ 
2 x 10 7 'A-^^^^m-^S'jy^t 3 / 4 . The maximum 
electron Lorentz factor is m e /m p lower and the syn- 
chrotron photon peak energy is w sy O^A^ 1 /? 2 2 keV, 
but due to the high Thomson depth these photons are 
thermalized to the blackbody value e 7jm ~ 0.1 keV. 

GRB protons and interactions in the SNR. — The GRB 
isotropic-equivalent 7-ray luminosity is L™° ss 10 52 L 52 
ergs/s. The corresponding GRB total proton luminos- 
ity is dN p /dt = L 1 lYm p c 2 = 2 x lO 52 !^^, s -1 . As- 
suming that a fraction £ p < 1 of these is accelerated 
in the internal shocks results in a proton distribution 
d 2 N p /dE p dt ~ 6 x 10 54 i 52 ^-Ep;G e v GeV-V 1 . High 
energy protons interact with synchrotron and inverse 
Compton scattered photons in the GRB fireball shock 
dominantly through pj — > A — > nir + /pn resulting 
in > 100 TeV neutrinos. The cross-section for A pro- 
duction at threshold E p E 1 ~ 0.2r 2 GeV 2 in the observer 
frame is ~ 0.1 mb. The optical depth for E p > 10 7 
GeV is t P7 > 1 in the high ~ MeV photon density in 
the internal shocks. Since protons lose ~ 20% of their 
energy per interaction, each proton will undergo a cou- 
ple of interactions, thus roughly half of the protons with 



E p > 10 16 eV will be converted to neutrons and escape. 
At lower energies, protons may be prevented by fireball 
magnetic field from escaping to the SNR. However, since 
internal shock radii range over 3 x 10 12 — 3 x 10 14 cm, 
this may allow a significant fraction of < 10 16 eV pro- 
tons to interact with the SNR shell. This implies that a 
fraction rj p (E p ) < 1 of protons escape the fireball shock 
region to propagate outwards. The isotropic-equivalent 
observer-frame proton luminosity impacting the shell is 

= 6 x W^E-^L^piEp) GeV-V 1 . (4) 

The large p 7 optical depths for both the SN-shock e sn 
and the pulsar MHD wind shock e m photons [Eqs. i ffl l 
and J2J] trapped inside the SNR shell causes most of the 
incoming protons above the thresholds [Eqs. (|TJ and J2J] 
to undergo photomeson production. Low energy protons 
below A production threshold may interact with SNR 
protons. If £ s h (typically < 10~ 4 ) is the fraction of shell 
protons which is accelerated to relativistic energies by the 
collisionless MHD shock going through the SNR shell, 
the column density of cold protons in the SNR shell is 
S p ~ I0 28 ( sh m snr vg 2 tj 2 cm -2 , where £ sh = 1 - 6h ~ 1- 

The total cross-section for pp interaction has been mea- 
sured up to very high energies (120 mb at y/s = 30 TeV 
^H) m accelerator experiments. We take the mean pp 
cross-section to be <cr pp >« 100 mb in the TeV and above 
energy range. The corresponding mean optical depth is 

<t pp >= S <cr pp >w l0 3 Cshm sn ,v 9 2 t^ 2 . (5) 

To calculate the neutrino flux from pp interactions, one 
needs to know the secondary charged particle multiplici- 
ties as a function of yfs. Analytic calculations of this mul- 
tiplicity, known as the KNO scaling law ^1 > are based on 
the rapidity (y = hx[(E+p z )/(E— p z )]/2, taking the beam 
direction along the z-axis) distribution of the secondary 
charged particles. Analytically the total charged particle 
multiplicity increases as ~ ln(s). Accelerator data show 
a slightly faster grow th in pp interactions measured up 
to y/s = 1.8 TeV]l3|. An extrapolation up to yfs ~ 400 
TeV, relevant for our calculation, for and mul- 
tiplicities from pp interactions using ~ ln(s) scaling is 
certainly very conservative. There exist in the literature 
other faster growing models such as ~ s 3 14]. The 
dominant neutrino production channels in our calcula- 
tion are pp — + ^ jK ± — ► fiv^ — > ev e v p v p . High energy 
photons are also produced from ir° — > 77. 

We used the PYTHIA 6.2 event generator 0] to sim- 
ulate pp interactions |l6|. The angular deviation of a 
secondary particle is related to its Lorentz invariant ra- 
pidity (y) as 9 ~ 1/ cosh(y). We select 7r ± and K ± which 
are forward, namely, y > 0. The average and 
multiplicities measured at yfs — 540 GeV in the pseudo- 
rapidity (77 ~ y for momentum p 3> m) region < r\ < 5 
are <n n >= 11.15 and <tik>= 1-25 [l3j. These numbers 
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agree with our PYTHIA simulations within error bars. 
For y > we have <n 7r >= 15.3 and <riK>= 1.7 at 
■s/s = 540 GeV from simulations. At the highest energy, 
E p = 10 20 eV (Vs = 4.33 x 10 5 GeV), ln(Vs) extrapo- 
lation gives <n n >= 15.3 x ln(4.33 x 10 5 /540) « 103, 
about 75% of our simulated value of 136. Similarly 
<n K >= 1.7 x ln(4.33 x 10 5 /540) » 11.4, also about 
75% of our simulated value of 15.2 at E p = 10 20 eV. 
For our calculation we used more conservative values: 
<n n >= 103 and <n K >= 11.4 at E p = 10 20 eV and 
la(y/s) scaling at lower energies. 

Secondary charged particles {ir and K^) in the pp 
interaction follow a X/E energy distribution (l6j . The 
energy of a particle of mass m ranges from m cosh(y)j cm , 
where 7 cm is the Lorentz boost of the center of mass in 
the lab frame and cosh(y) is the rapidity dependent boost 
factor, up to the primary proton energy E p . For y > 
we use for secondary and the range 



7cm"V,lC 



(6) 



The pions and kaons decay into neutrinos with Lorentz- 
expanded decay times {t.k,k1-k,k) in the observer frame. 

In the tt^ / — > eUeU^V^ decay the e and 3 z/s each 
share roughly 1/4 of the pion/kaon energy. The branch- 
ing ratio for the kaon decay channel is ~ 64%. The v 
multiplicity for an incident proton of energy E p < 10 11 
GeV, using Eq. JJjJ, is then 
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GeV ~ 4 GeV 



(7) 



for each type of neutrinos: v e , and v^. Here is a 
step function following from Eq. (j^J. The normalization 
factors N-„ and Mk are found by integrating the 1/E dis- 
tribution of tt^~ and K ± in the energy range given in Eq. 
© for E™\ sa 10 11 GeV and equating to the respec- 
tive total numbers, J\f n =<n 7r > /In (£'™ x /TO 7r 7 cm ) w 7, 
M K = 0.64 <n K > /In (E^/m Klcm ) a 0.6. 

An additional pp component may result from protons 
accelerated by the collisionless shock in the shell col- 
liding with shell thermal nucleons. For a typical SNR 
shell energy lO 51 ^! ergs this neutrino component has 
< ^s/jlO 51 ^! ergs, at least an order of magnitude below 
the GRB shock proton component, but it could become 
important for exceptionally energetic supernova shells. 

Neutrino flux calculation. — The GRB protons from 
internal shocks undergo pp interactions below E p t h from 
Eqs. (JIJ and (J2J, and pr) interactions above E p ^h- The 
neutrino flux = d 2 N /dE v dt) at Earth from a single 
GRB-SNR at distance D is 



$>u = 



1 / / f P p M v{E p )-^- t dE p \ E p <E p ^ th 



4^D 2 1 (/„/4) 



d 2 N 
dE p dt ' 



for each v type. Here f pp = min(l, <T PP >) from Eq. 
(JSJ, and ./vis the fraction of E p lost to 7r's in pr) in- 
teractions |2j, which is ~ 1 for an SNR with r pl 1. 
The synchrotron cooling energy E s v is found by equating 
the decay times (tj7j with j — n, K, ji) and synchrotron 
cooling times for 7r's, K's and /i's. For the magnetic field 
in Eq. we get the maximum energies as 

» 2.3 x 10S»JVVWOT 0) 

where rn^'s are in GeV and Tj's are in s. Taking m snr ~ 
10 and typical age t^ ~ 50 d, the maximum energies 
from Eq. @ are E* ~ 10 10 GeV, E S K ~ 10 11 GeV and 
El ~ 10 9 GeV, where we used £b w 10~ 2 £_2 as in typ- 
ical GRB fits. The corresponding v steepening break 
energies are > 10 9 GeV. The cooling times for inverse 
Compton scattering are longer than respective lifetimes 
in the Klein-Nishina limit for td > 10 d. 

The detection probability for a muon neutrino in an ice 
detector is V « 1.7 x lO" 6 ^ /TeV) 13 , where [3 = 0.8 for 
E v < 1 PeV and /3 = 0.36 for E v > 1 PeV. Multiplying 
this with the flux of Eq. (JSJ and the burst duration, a 
numerical integration gives the expected event rates in 
a km 2 detector from a single GRB. The total number 
of ^-induced upward muon events for a single GRB of 
luminosity 10 52 L52 ergs/s and burst duration At = 10Ai 
s at D r> 1.6 x lO 28 ^ 1 ^ + z)/2 - (1 + z) 1 ^ 1 / 2 ] cm 
~ 10 28 2 Z?28.2 cm for redshift z ~ 1 is 



V 
V 



TeV-PeV 



PcV-EoV 



= l.l£)^ 8 2 2 Aie P ?7 P (S)£52 

= 0.06^ 28 2 2 A 1 ^?7 P (S)^52. (10) 



E p > E Pi th 



Here we assumed f pp — 1 and f n ~ 1 for m snr ~ 10 and 
10 ^ 100 days. These single-burst numbers can be 
~ 10 2 times larger, e.g. for rare bursts with D = 10 27 5 
cm and At = 100 s occurring at rates ~ 3 yr _1 . 

The diffuse i/-fhrx is obtained from Eq. JSJ multiplied 
by At and an observed rate ~ 600/47T yr^ 1 sr _1 . Figure 
^shows the diffuse v p energy flux per decade (same for v e 
and Dfj) from GRB models with pre-ejected SNR shells. 
We assumed ~ 10 _1 of protons escape the fireball shocks 
to reach the shell, where f pp = 1 and Ef, = 10 9 GeV. 
The top curve is calculated assuming all the GRBs have 
a SNR shell. This is an upper limit for maximum j/'s 
through pp and pj interactions in the SNR shell. The 
bottom curve assumes 10% of the GRBs have an SNR 
shell. Also shown are the diffuse ^-fluxes from 07 inter- 
actions in the internal shocks of bursts (WB [2|, short 
dashed curve), and in GRB afterglows (|3j, short dashed 
curve), as well as the ^-flux obtained from the cosmic ray 
limits (WB Limit i° n g dashed straight lines). In a 
km 2 detector 0], the number of diffuse events assum- 
ing 10% of GRBs have SNR shells is ~ 6 yr" 1 sr" 1 at 
TeV-PeV, and - 0.3 yr" 1 sr" 1 in the PeV-EeV range. 

Discussion and Implications. — The single-burst in- 
fluxes calculated here are predicated on the existence of 
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same GRB-contemporaneous internal shock protons in- 
teracting with thermal 0.1 keV photons within the shell 
wall, whereas consider such protons interacting with 
photons from the MHD wind inside the shell cavity. 

The pre-ejected supernova (supranova) model of GRB 
is a subject of interest and debate yjj for interpreting 
the 7- and X-ray data, and independent tests would be 
useful. The neutrino fluxes discussed here provide such 
a test, the predicted event rates being detectable with 
kilometer scale planned Cherenkov detectors. 

SR thanks J. P. Ralston for helpful discussions. Work 
supported by NSF AST0098416. 



FIG. 1: Diffuse neutrino flux (Ef,$v) from post-supernova 
(supranova) models of GRBs (solid curves), assuming that 
(top curve) all GRBs have an SNR shell, or (bottom) 10% 
of all GRBs have an SNR shell, and 10 _1 of the fireball pro- 
tons reach the shells. Long dashed lines correspond to the 
Waxman-Bahcall cosmic-ray limit, short dashed curves are 
the diffuse v flux from GRB internal shocks and afterglows. 



a pre-ejected SNR shell by the progenitor of the GRB, 
which occurs at the same location after a delay of weeks 
||. The j/s produced by pp and p7 interactions between 
GRB relativistic protons and SNR shell target protons 
and photons will be contemporary and of similar duration 
as the GRB electromagnetic event. The high pp optical 
depth of the shell also implies a moderately high average 
Thomson optical depth tt oc t^ 2 of the shell, dropping 
below unity after ~ 100 d. Large scale anisotropy as 
well as dumpiness of the shell will result in a mixture of 
higher and lower optical depth regions being observable 
simultaneously, as required in the supranova interpreta- 
tion of X-ray lines and photon continua in some GRB 
afterglows [4j. Depending on the fraction of GRB with 
SNR shells, the contributions of these to the GRB diffuse 
zz-flux has a pp component which is relatively stronger at 
TeV-PeV energies than the internal shock pr) component 
of 0, and a shell pj component which is a factor 1 (0.1) 
of the internal shock pj component (Fig. ^) for a frac- 
tion 1 (0.1) of GRB with SNR shells. Due to a higher 
synchrotron cooling break in the shell, at E v > 10 17 eV 
the shell component could compete with the internal Q 
and afterglow 3| components. 

Our pp component is caused by internal shock- 
accelerated power-law protons contemporaneous with the 
GRB event, differing from Q who considered quasi- 
monoenergetic 7 P ~ 10 4 5 protons from an MHD wind 
over 47T leading to a~ 10 TeV neutrino months-long pre- 
cursor of the GRB. Our pj component arises from the 
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